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We consider the quantum dynamics of a harmonic oscillator in noncommutative space under the 
influence of linearized gravitational waves (GW) in the long wave-length and low- velocity limit. 
Following the prescription in^ we quantize the system. The Hamiltonian of the system is solved by 
using standard algebraic iterative methods. The solution shows signatures of the coordinate non- 
commutativity via alterations in the oscillation frequency of the harmonic oscillator system from 
its commutative counterpart. Moreover, it is found that the response of the harmonic oscillator 
to periodic GW, when their frequencies match, will oscillate with a time scale imposed by the NC 
parameter. We expect this noncommutative signature to show up as some noise source in the GW 
detection experiments since the recent phenomenological upper-bounds set on spatial noncommu- 
tative parameter implies a length-scale comparable to the length-variations due to the passage of 
gravitational waves, detectable in the present day GW detectors. 

PACS numbers: ll.lO.Nx, 03.65.Ta, ll.lO.Ef, 04.30.Nk, 42.50.Dv 



Gravitational waves (GW(s)) are tiny vibrations in 
spacetime structure itself. The present day scenario of 
GW detection experiments primarily consists of ground- 
based (LIGO\ VIRGO^, GE03, TAMA" etc.) and 
space-based (LISA^) interferometers. The key idea here 
is to measure the relative optical phase shift between the 
light paths in tvifo perpendicular km-length arm cavities 
where this phase shift is due to the relative displacement, 
induced by a passing GW, of the two mirrors hung at 
the end of each arm cavity. However, the history of ex- 
perimental GW physics began with resonant-mass detec- 
tors, pioneered by Weber in the 60's. In the following 
decades, although the sensitivity of resonant-mass detec- 
tors have improved considerably, it is clear that it could 
only allow the detection of relatively strong GW signals 
in our Galaxy or in the immediate galactic neighbour- 
hood. Nevertheless, building the interferometric detec- 
tors takes many years of preparation and considerable 
finance whereas the resonant detectors, being relatively 
small-scale instruments, are more easily realizable. Beis- 
des, the study of resonant- mass detectors is instructive in 
itself because it focuses on how a GW interacts with an 
elastic matter causing vibrations with amplitudes many 
order smaller than the size of a nucleus. In a resonant bar 
it is possible to measure vibrations which corresponds to 
just a few tens of phonons^, and variations AL of their 
length L, with -^ ~ 10~^^. It is, therefore, at the quan- 
tum mechanical level, that experimental evidence for the 
GWs is likely to appear—. 

Interestingly, in recent developments of noncommuta- 
tive (NC) quantum mechanics^"— and NC quantum field 



theoryi^— , where the coordinates x^ satisfy the NC al- 
gebra 



ie"" 



(1) 



the upperbounds on various NC parameters appear- 
ing in the literature^"— are quite close to this length 
scale. A wide range of theories have been constructed 
in a NC framework including various gauge theoriesi^, 
gravity^ and even encompassing certain possible phe- 
nomenological consequences^^"— . The upperbound on 
the value of the coordinate commutator 9"^^ found inii 
is ^ (lOTeV) which corresponds to 4 x 10~'*°m^ for 
h=c=l. Whereas such upperbounds on time-space NC 
parameter 6*"* is < 9.51 x 10~^^m^. However, recent stud- 
ies in NC quantum mechanics revealed that the NC pa- 
rameter associated with different particles are not same^ 
and this bound could be as high as 9 ^ (4GeV)~ — 
These upperbounds correspond to the 



(30MeV)~^^. 
length scale '-- 



10- 



-20^ 



10 



-17^ 



With the prospect of the direct detection of GW(s) 
of such tiny amplitude as '^ 10~^^ in the near future, a 
good possibility of detecting the NC structure of space- 
time would be in the GW detection experiments as it 
may as well detect the signature of noncommutativity. 
As a first step towards this endeavour, we have analysed 
the interplay of classical GW(s) with a free test particle 
in a NC quantum mechanical framework in"^ . Our analy- 
sis suggested that investigating a NC harmonic oscillator 
system interacting with the passing GW is more likely to 
reveal NC effect comparable with the effect of GW. More- 
over, a simplistic consideration of the interaction of GW 



with the resonant bar shows that the fundamental mode 
of vibration of a thin cyhndrical bar is formahy identical 
to a harmonic oscillator driven by a force exerted by the 
GW~. Therefore, NC quantum mechanical consideration 
of the interaction of GW with matter, specifically the 
harmonic oscillator which is essentially inherent in the 
resonant mass detector system, would certainly prove in- 
structive. In light of all these facts, we, in the present 
paper, formulate the quantum mechanics of a NC har- 
monic oscillator, interacting with a linearised GW in the 
long wave-length and low velocity limit. 

Since it has been demonstrated in various formulations 
of NC general relativitySi,^!,^^ that any NC correction in 
the gravity sector is second order in the NC parameter, 
therefore, in a first order theory in NC space, the lin- 
earised GW remains unaltered by NC effects and any NC 
correction appearing in the system will be through the 
matter part only. This is true not only with the canoni- 
cal (i.e. constant) form of noncommutativity but also for 
the Lie- algebraic NC spacetime'^^. We shall incorporate 
the NC effect by writing the NC Hamiltonian for the sys- 
tem and then reexpressing it in terms of the commutative 
coordinates and their momenta by the well known Bopp- 
shift transformations'^^. This commutative equivalent 
model will be quantized following—. We shall deal with 
the linearly polarized GWs in this analysis for simplicity. 

In a linearized theory of gravity the connection and 
curvature tensor take the form^ 

r^'a/S = -77^" ihfi^,a + ha„,i3 - h^p^u) (2) 

^afSfiv ^^ 7T v'^az^,/t/3 "r f^Pfi.au '^Pu.a^ f^afi-Pu) ['^J 

where /i^,^ is the metric perturbation on the flat 
Minkowski background 77^1,. Next we choose the 
transverse-traceless (TT) gauge 



is its mass and the GW is treated as an external classical 
field. Note that, eq.® can be used to describe the evo- 
lution of proper distance in TT-gauge frame as long as 
the spacial velocities involved are non-relativistic. Also, 
\x^\ has to be much smaller than the typical length scale 
over which the gravitational field changes substantially, 
i.e. the reduced wavelength -^ of GW. The above condi- 
tions are collectively referred to as the small-velocity and 
long wavelength limit. Thus, with eq.® we can analyse 
the interaction of GW with a detector which has a char- 
acteristic linear size L <^ ^. Note that this condition is 
satisfied by resonant bar detectors as well as earth bound 
interferometers but not by the proposed space-borne in- 
terferometers such as LISA^ or by the Dopplcr tracking 
of spacecraft. 

The Lagrangian for the system, whose time evolution is 
described by eq. ([6]) , can be written, upto a total deriva- 
tive term'^^ as 



L — —mx 
2 



mT^okijX 



22 
—mvj X 
2 



(7) 



The canonical momentum corresponding to Xj is pj 
mxj — TTiFQj.a;'^ and hence the Hamitonian becomes 



''=2^y^^ 



mVl^x 



'- 22 
—mzu Xj 

2 ^ 



(8) 



Assuming that the GW is propagating along the 2;-axis, 
due to the transverse nature of GW(s), F-'ofc has non-zero 
components only in the x — y plane and therefore the 
response of the pendulum to it is essentially confined to 
that plane. To impose noncommutativity on this plane, 
we assume that the coordinates follow the algebra ([1]) and 
"quantise" the system on this NC plane. To proceed, we 
replace x^ and Pj in the above Hamiltonian by operators 
x^ and Pj satisfying the NC Heisenberg algebra^ 



/iom = , 



V;^ = 



K = ^ 



(4) 



[XuPj\ 



ihSi. 



[xi,x.j\ = iOe^j , [pi.pj] = . (9) 



to remove all the gauge redundancies of the theory 
and the GW is characterised by the only non-zero 
components^ hn = — /122 and /112 = /121, called the -f and 
X polarisation respectively. The only non-trivial compo- 
nents of the curvature tensor in TT-gauge are^ 



R''o,ko = 



dt 



l(fhjk_ 
2 dt^ 



(5) 



and the geodesic deviation equation in the proper detec- 
tor frame becomes^ 



d^x^ 



-m.R^o^kox'' 



mzu X-' 



(6) 



Here t is the coordinate time of the proper detector frame 
and is same as it's proper time since we are confining our- 
selves to first order in the metric perturbation. Eq.® 
governs the response of a 2-dimensional harmonic oscil- 
lator with frequency w to the passage of a GW. Here x^ 
is the proper distance of the pendulum from the origin, m 



It is well known that this can be mapped to the stan- 
dard {0 — 0) Heisenberg algebra spanned by Xi and Pj 
using2ii2^ 



X,, 



^^^..^., P^ 



P, 



(10) 



Using the traceless property of the GW and rewriting the 
NC version of eq. ^ in terms of the operators Xi and Pj , 
we obtain 



H 



P,,^ 1,9, mtu^ 

2to 2 J ^ ofc J /^ 2h 



]m^ 



2% 



f ■ P P,V^ 



(11) 



In the above equation the first two terms are for the or- 
dinary harmonic oscillator, the third term, linear in the 
affine connections, shows the effect of the passing GW on 
the ordinary harmonic oscillator system, the fourth term 
is the signature of NC space, a pure NC term linear in 



the NC parameter and the final term shows the couphng 
between the GW and spatial noncomniutativitj*^. Since 
the length-scale range describing various upper-bounds 
on NC parameter that we have discussed earlier is quiet 
close to the detectable variations of length AL in the 
GW detection experiments, at least the fourth term in 
the left hand side of equation ([TT|) should be compara- 
ble with the third (purely gravitational) term and show 
up in the GW experiments as some noise source. This, 
in essence, is similar to some recent results^^ where the 
coordinate noncommutativity is shown to generate holo- 
graphic noise claimed to be detectable in the signals of 
the interferometric GW detectors GEO6002£. 

Defining raising and lowering operators in terms of the 
oscillator frequency zu 



X, 



Pj = 



f ^ 



V 2mw 



1/2 



f hmvj\ 



1/2 



(12) 
(13) 



we write the Hamiltonian pT|) as 
H — fiTu I ajttj + 1 



ifi ■ 



J J 



8 



— -mw 



■dejka]jak 



(14) 



where C.C means complex conjugate. Working in the 
Heisenberg representation, the time evolution of aj (t) is 
given by 



dt 



-iwa 



] 



1 ■ 4. mvj'^9 

2 -^ " 2h 
imzud / ; . ; \ / f ' 



-^jkO-k 



H — ^r— [dijhik + eikhijj I flfe - aM (15) 

and that of a'At) is the C.C of the above equation. Next, 
noting that the raising and lowering operators must sat- 
isfy the commutation relations 



[aj{t),ak{t)] =0= \a]{t),al{t) 



(16) 



we write them in terms of aj{0) and a-{0), the free op- 
erators at time t = 0, by the time-dependent Bogoliubov 
transformations 



u^v , uu' — vv'' = I, written in matrix form where T de- 
notes transpose, j denotes complex conjugate transpose 
and / is the identity matrix. Since aj{t = 0) = aj(0), 
Ujk(t) and Vjk(t) have the boundary conditions 



,fc(0) = / , v,k{0) = 



(18) 



Then, from eq. (|15p and its C.C, we get the following 
equations of motions in terms of (^ = u— w^ and f = u+v^: 

-^ = -iTJjijk - ^hjiCik ^^^i'C;fc (19) 

d^jk ^ 1 ; ^ ^ ^ mzu'^9 

-^ = -i^Qjk + -;^hji£,ik + QjiCik ^^^J'^'fe (2^) 

where Qji is the term reflecting the interplay of noncom- 
mutativity with GW 



e, 



imw9 
Ah 



\ lljm^jnl ^jni'^ml 1 ■ l^^j 



Eq(s) p9l20p are difBcult to solve analytically for general 
hjk- However, our goal, in the present paper, is to inves- 
tigate whether we get comparable effects of spatial non- 
commutativity and gravitational wave on the harmonic 
oscillator system in the simplest of settings. Therefore 
we shall solve eq(s) (|19l20p for the special case of lin- 
early polarized GW(s). 

In the two-dimensional plane, the GW, which is a 2 x 2 
matrix hjk, is most conveniently written in terms of the 
Pauli spin matrices as 

h,k (t) = 2 fit) (exfTJ, + e+a;\) = 2f{t)sA<jf, . (22) 

Note that the index A runs from 1 — 3, however, no contri- 
bution from cr^ is included. 2f{t) is the amplitude of the 
GW whereas Ex (t) and e+ (t) represent the two possible 
polarization states of the GW and satisfy the condition 
el- +e'^ = 1 for all t. In case of linearly polarized GW(s) 
however, the polarization states sa are independent of 
time and f(t) is arbitrary. To set a suitable boundary 
condition we shall assume that the GW hits the particle 
at i = so that 



fit) = , for i < 0. 



(23) 



We now move on to solve eq(s) (|19I20I) by noting that 
any 2x2 complex matrix can be written as a linear com- 
bination of the Pauli spin matrices and identity matrix. 
Hence we make the ansatz : 



Ojit) = Ujkit)akiO) + Vjkit)aliO) 
a] it) = aliO)ukjit)+akiO)vkjit) 



(17) 



Qk it) 

ilk it) 



AI,k 

Cljk 



Di<y]k 



B^a' 



D2<J^ 



jk- 

.2 
jk 



B,a% 



Dsa 



]k 



(24) 
(25) 



where the bar denotes the C.C and Ujk and Vjk are the 
generalised Bogoliubov coefficients. They are 2x2 com- 
plex matrices which, due to eq. (I16p . must satisfy uv'^ = 



Substituting for hjk, (jk and ^jk from eq(s) ([^ . ([M)) 
and (psj) in eq(s) (|19l20p and comparing the coefficients 
of / and (T-matrices, we get a set of first order differential 



equations for A, _Bi, i?2, ^3, C, Di,D2,D3 : 



where 



A = -ivjC - f {eiBi + 6383) - iAB2 

Bi = -ivoDi- f{eiA-ie3B2)+KB3 

B2 = -iiuD2 - if i^sBi - eiBs) - iAA 

B3 = -imD3-f{e3A + ieiB2)-ABi 

C = -iwA + f{eiDi+e3D3) + 4iXf{e3Bi-eiB3) 

-iAD2 
D, = -iwB, + f{eiC- 16302) +4Xf{ie3A^eiB2) 

+AD3 
D2 = -imB2 + ii{e3D^~eiD3) + A\j{eiB^+e3B3) 

-iAC 
D3 = -izuB3 + f{e3C + ieiD2)-4\f{ieiA + e3B2) 

-ADi (26) 



where A and A are given by 



A 



m,-cu 

2h 



A 



Ah 



(27) 



and the dot represents derivative with respect to time 
t. Noting that \f{t)\ « 1, the above set of equations 
can be solved iteratively about its f{t) =0 solution by 
applying the appropriate boundary conditions P^ [^ . 
We therefore obtain to first order in the gravitational 
wave amplitude and also in the NC parameter 



A{t) = C{t) = e-*^* cos(Ai) 



B2{t) = D2{t) = -ie 
Bi{t) = -D,{t) = -e 



■ sin(Ai) 



(28) 
(29) 



[(eiCos(At)-e3sin(Ai))/(i) 



2injei / di'e*"(*-*') cos(At')/(i') 



-2iwe-j 



di'e^^(*-*') sin(Ai')/(t') 



+w' 



(ei f gi{t')dt' -83 J g2{t')dA (30) 



B3{t) = -D3{t) = -e 



[(£3 cos(At) + ei sin(Ai))/(i) 



+2iws3 / di'e*^(*-*') cos(At')/(i') 

JQ 

+2iwei [ dt'e'"^^'-''^ sm{At')f{t') 
Jo 

+^2 (£3 I g,{t')dt' + £1 /" g2it')dA (31) 



5i(i)= / dt'e-'^'' cos{At')f{t') 
Jo 

g2{t)= f dt'e-'^*' sm{At')f{t') . 
Jo 



(32) 



The system has now been essentially solved once we spec- 
ify the initial expectation values of the pendulum's posi- 
tion rg = (a;o,yo) a-nd momentum p^ — {pxa,Pya) when 
the GW just hits the system at time t = 0. The time 
evolution of the coordinates can be calculated employ- 
ing the following scheme. Combining the expressions for 
A,Bi,B2,B3,C,Di,D2,D3, we can write the solutions 
for C and ^ using Eq.(s) (|24l25p which in turn give u 
and V. Using Eg. dlTp . we can now combine u and v into 
aj{t) and ajt). From the initial position and momen- 
tum expectation values, i.e.(ro) — (Xi(0), ^"2(0)) and 
(Pq) — {Pi{0),P2{0)), we get the raising and lowering 
operator qj (0) and a- (0) at time t = 0. We then use 

them in Eq.(s) p7)) to find aj (t) and a, (t) at a general 
time t and these yield the time evolution of the expecta- 
tion values of position coordinates {Xi (t)) and {X2 (t)) 
of the pendulum. The general expression of {Xi (t)) thus 
obtained is given by 



{X,{t)) = [Re{A)+Re{D*)]X,{0) 
+ [Im(D2) + Re{D\)]X2 (0) 
,A(0) 



+ [-Im{A) + Im{Dl)] 
+ [Re{D2) + Im{Dl) 



m-cu 

P2{0) 



(33) 



where Re and Im denote the real and imaginary parts of 
the A, Di, D2 and D3. Substituting their values from the 
equations ^, ^, ^ and ^ will give the expres- 
sion for {Xi (t)) for a general GW amplitude f{t). Since 
it will be difficult to see through this much complicated 
expressions we assume a monochromatic GW waveform 
oscillating with frequency m' , f (t) — /oe'^ * which, upon 
substitution in equation f (P51) gives 



(^1 (t)) 



(cosn7_i + costu+i) , . {smzu^t + sinvj+t) (cos w+t — cos -cu^t) 
n ^i(O) + TT-- n(0) + A2(0) 



+ 1 



2mw 



(sin ra7_|_t — sin w-t) 
2mw 



/o{eiXi(0) + £3^2(0)} 



(ei + £3) /o 



+ (ei - £3) /o 



Acc7_ 



An7_t / , , Am^t Pi(0) Atn^i 

sm ^-— Xi 0) sm ^ cos ^-— 

2 V 2 TnciJ 2 



zu'_ V Auj.t ( , , AvoA P2(Q) Aw+t 
sin ^^ ( X2 (0) sin -—^ ^^ cos ^^ 



A 



VJa 



mw 



tuL \" . Atn+t / ^ ,„, . Atn+t Pi(0) Aw+t 
sm — ; — I Ai(0) sm — cos — - — 



Aro 






Aot- 



2 

AvD^t 



sm ■ 



. An7„t P2(0) Aw.t 
X2 (0) sm — cos — - — 



mw 



where Aw = zu — tu' . zu± — tu ± A, w'^. = n?' ± A and Atu± — Aw ± A. 

I 



P2(0) 



(34) 



This result impHes that the presence of noncommutativ- 
ity alters the response of the harmonic oscillator to a 
periodic GW from its commutative counterpart. When 
the frequency of the GW is very close to that of the har- 
monic oscillator (Aro ~ 0), the oscillatory terms present 
in the solution representing the response of the system 
to the GW will oscillate with frequency A with a large 
amplitude. It should be possible to detect this effect. 
Now putting = 0, i.e. A = gives us the classical 
result of GW interacting with a harmonic oscillator in 
the low-velocity, long-wavelength limit whereas putting 
/o = 0, i.e. in the absence of gravitational wave the 
solution assumes the form of a NC harmonic oscillator. 
Interestingly, the presence of |- factor in A, i.e. in the 
NG correction terms even after the computation of the 
expectation value indicates that the NG effect is inher- 
ently quantum mechanical in nature. Similar expression 
for (X2 (t)) can also be obtained. Further realistic sce- 
narios can be obtained if we use various forms of periodic 



GW signals with more than one frequency and do similar 
computations. Now that we have studied the interaction 
of a single NC harmonic oscillator with GW, a further 
advancement will be to extend it for a macroscopic piece 
of elastic matter. In fact in a resonant mass detector it is 
possible to detect vibrations which are incredibly small, 
with amplitude many orders of magnitude smaller than 
the size of a nucleus. In that context studying the ef- 
fect of noncommutativity may prove interesting. Work 
in this direction is in progress and will be taken up in 
the subsequent papers. 
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